Aharonov-Bohm Oscillation and Chirality Effect in Optical Activity of Single Wall 

Carbon Nanotubes 
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We study the Aharonov-Bohm effect in the optical phenomena of single wall carbon nanotubes 
(SWCN) and also their chirality dependence. Specially, we consider the natural optical activity as a 
proper observable and derive it's general expression based on a comprehensive symmetry analysis, 
which reveals the interplay between the enclosed magnetic flux and the tubule chirality for arbitrary 
chiral SWCN. A quantitative result for this optical property is given by a gauge invariant tight- 
binding approximation calculation to stimulate experimental measurements. 
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Aharonov-Bohm (AB) effect manifests the significance 
of the global nature of the vector potential in quantum 
theory 1,2 . Such a geometric phase effect will result in a 
conductance oscillation with respect to the enclosed flux 
in the cylindrical conductor. This phenomenon is known 
as the Altshuler-Aronov-Spivak effect. Since the dis- 
covery of carbon nanotube 4 , it provides an ideal hollow 
cylindrical lattice sheet with distinguished chiral struc- 
tures which stimulates extensive studies*** 7 **^ in recent 
years. Moreover it provides a novel system for studying 
the interplay between the enclosed AB magnetic flux ver- 
sus the chiral symmetry of the tubule. When an applied 
magnetic flux is threaded through, it has been shown el- 
egantly in the tight binding approximation (TBA) (or 
effective mass approximation)***-!*!^ that the fundamen- 
tal gap of a single wall carbon nanotube (SWCN) is a 
periodic function of magnetic flux. This result is also 
referred to a kind of AB effect. The corresponding spec- 
trum feature has been further analyzed in Ref*-! where 
the periodicity of Van Hove singularities with respect to 
the magnetic flux was also addressed**. Experimental 
transport studies* 3 ** 4 *** at the finite temperature, for 
multi-layer carbon nanotubes, also show consistently a 
current oscillation with a single flux quantum </)q . 

Since the SWCN may behave as either an insulator or 
a metallic conductor depending on its chirality as well 
as the strength of penetrating magnetic flux, the con- 
ventional transport study might not be appropriate for 
exploring the generic chirality dependence of the AB ef- 
fect specially for the insulating cases. Actually, as long as 
the AB flux appears in the hollow SWCN, all electronic 
canonical momenta will acquire a corresponding vector 
potential. As a result, the exponential phase factor as- 
sociated with each C-C link induced by the threading 
magnetic flux together with the orientation of the chiral 
trident will make for an interplay between the threading 
flux and chirality. Therefore the AB-type effect should 
appear in a variety of phenomena such as the optical 
properties other than the transport measurements. 

In this paper, we show the AB oscillation in the natural 



optical activity (or natural gyrotropy) of the SWCN with 
arbitrary chirality. We derive a generic expression for 
the natural gyrotropy based upon a systematic symmetry 
analysis and further complement a gauge invariant TBA 
calculation. It exhibits the interplay among the light 
polarization, tubule chirality and threaded magnetic flux 
and reveals the characteristic role of the chiral index 17 . 




FIG. 1: For given Ro[ni,n2], there are five chiral vec- 
tors, Ri[m + 7i2,— m], R'2[n2,~ m - n 2 ], Ra[— hj.,— ria], 
R4[— m — ri2, Tii] and Rs[— 712, ni + 712], denoting exactly the 
same SWCN with their chiral angle deferred 7r/3 from each 
other successively. The tubule described by Rq[tii +712, —712] 
is the mirror image of that described by Ro. ai and a2 are 
2D graphite lattice basis vectors. 

An ideal SWCN can be viewed as a graphite sheet 
rolled up into a tubule in a fixed way along the chi- 
ral vector R = niai + 71282 with ai and a.2 being the 
2D graphite lattice basis vectors. Its geometric struc- 
ture is hence described by a pair of integers [ni, 712]. For 
given [ni, 712], there are six chiral vectors, i.e., R,o[m, ^2], 
Ri[ni + ri2, — n±], Rs[— ?i2,ni + n^] as shown in 
Fig. 1, denoting exactly the same SWCN with chiral an- 
gles (the angle between R and ai) as 0, 8 + 7r/3, . . ., 
6 + 57r/3 respectively. If we introduce the chiral index 
v = mod[?ii — 7i2,3] as in Ref. 17 , we find that they can 
be divided into two subsets, one is {Ro, R2, R4} with the 
same v and the other is {Ri,R3,Rs} with an opposite 
v. The essence of this obvious but nontrivial property, 
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the same SWCN with opposite chiral indices, lies in the 
sameness of the two in-cell carbon atoms A and B. Ac- 
tually this is the heritage to the carbon tubule inherited 
from the 6mm symmetry of the graphite sheet which has 
the 3m symmetry as its invariant subgroup with corre- 
sponding quotient groups as E and C2(or a). Ci (or a) 
will reverse the A and B atoms and also the sign of chi- 
ral index. As a direct result of the above analysis, the 
physical quantities of SWCN can be written as a peri- 
odical function of 9 with period 27r/3 for v — ±1 and 
7r/3 for v = 0, which are the consequences of the 3m 
symmetry and 6mm symmetry, respectively. When per- 
forming Fourier expansion for these physical quantities 
with respect to 9, all the expansion coefficients should be 
functions of those invariant qu antities such as t he length 
of chiral vector A|ai| (A = yj n\ + n\ + nin 2 ) and the 
greatest common divisor of n\ and ri2 (denoted by N). 

Since the rotation of the linear polarized light trav- 
elling through SWCN along the tubule axis carries the 
information of its chiral structure and the enclosed mag- 
netic flux, we consider the optical rotation power(ORP), 
i.e., the rotation angle per unit length, as an observable to 
investigate the interplay between chirality and threaded 
magnetic flux for the SWCN. By fixing the direction of 
the incident light with frequency u parallel to the tubule 
axis (z-axis), the ORP has the expression 



X = u 2 g xyz /(2c 2 ) 



(1) 



where c is the light speed and the third rank tensor g xyz 
is the derivative of the xy component of dielectric tensor 
with respected to q z in the long wavelength limit q —> 0, 
i-e., g xy z = de xy /(idq z )\ q ^(r&. Notice that a SWCN and 
its mirror image with respect to the horizontal plane pos- 
sess the same chiral index, but different chiral angle 9 and 
—9, respectively (see Fig. 1). Further taking considera- 
tion of the third rank tensor properties, one can see that 
\ v should be an odd function of 9. Hence, it has a Fourier 
series with only odd parity sine functions 



x ± (9) = af sin(36>) + af sin(66>) + af sin(90) + 
X°(9) = a£sin(6(9)+a2sin(1260 + -- - , 



(2) 



where the superscripts ± or denote the chiral index, and 
the coefficients a v n depend on the curvature, the magnetic 
flux and the frequency of incident light. As a natural 
consequence of Eq. {2J the ORP of zigzag (9 — 0) and 
armchair (9 = 7r/6) tubules is zero. Moreover, consider- 
ing the rotation of R by 7r/3, since Ro and Ri correspond 
to the same SWCN with opposite v, it is straightforward 
to see x ± {Q + t/3) = X T (#)- Hence, we have 



o*(w,A;^) = (-ir^(o;,A 



(3) 



We further adopt the scheme developed by White et 
ali&, which exhibits the chiral structure of SWCN ex- 
plicitly with cylindrical geometry being properly built in. 
In this scheme, all the lattice sites can be generated by 
repeating the pure rotation Cjv and the screw operation 



S(a,h), while the latter one is shifting along the tubule 
axis by h with a simultaneous rotation around the tubule 
axis by a. h and a can be obtained through equation 
h = V3N\a 1 \/2A, a = (2p 1 n 1 + 2p 2 n2+pin 2 +P2n 1 )TT/A 2 
with integers pi,p 2 satisfying p 2 ni — pini = N . Hence, 
the Bloch momenta n £ [— tt, 7t) and n = 0, 1, • • • , N— 1 as 
good quantum numbers can be extracted from the char- 
acters of the U (l)-representations of S(a, h) and Cjv, re- 
spectively. If the magnetic flux is threaded through the 
SWCN, the space displacement groups S(a, h) and Cjv 
will be replaced by the corresponding magnetic displace- 
ment groups, which are again commutable to the Hamil- 
tonian of the flux threaded SWCN. The Bloch momenta 
then become n + 4>/4>o and k + a{<j)/4>o), respectively. 
This provides the generic way of the flux dependence for 
all the gauge invariant energy spectra^ as well as ma- 
trix elements of physical quantity. As a result, the flux 
dependence of the observable for intrinsic SWCNs will 
exhibit an AB oscillation with a flux period tf>o. This can 
be easily seen from the fact that, these quantities can be 
often expressed as a double summations over k and n, 
the former sums over [— tt, tt) and then smears the flux 
dependence via k + a(<p/<po) while the latter summation 
will contribute the AB type flux dependence. Then we 
conclude that a v n in Eq. (2) are periodic functions of <j> 



<(^)A; 4>) = fl£(cJ,A;</> + 0o) 



(4) 



Equations (J2J , © and complete the generic symmetry 
analysis for the OPR. 

Now we apply the TBA to the SWCN with arbitrary 
given pair of [711,712], and ignore the curvature effect for 
convenience. The TBA Hamiltonian for the SWCN with 
threading flux along the tubule axis can be written as a 
2x2 matrix form in the momentum space 



H(k, n) = Vq 





7n(«) 



(5) 



with Vo being the transfer integral equal to 2.6 eV. The 
off-diagonal matrix element has the form 7 n («0 = 1 + 
e<* with/?! = ^ +a ±y^{ n+ ±) and/^ = 
■). Then it is straightforward 



to obtain the energy spectrum En'' v \n) = ±Vb|7n(ft)|^ 
and the one-particle wave function for the conducting or 
valence band reads 



\.,nAv)) = ±{\^n,A)±^ 

V2 V l7n(«)| 



K, 71, B\ 



(6) 



where the plus and minus signs correspond to the con- 
ducting band and the valence band, respectively. The 
Bloch sum in Eq. © reads 



\k,ti, s) 



1 N 
V2MN f-f 

I — 1 



M 

E 

■i=-M 



e lKm+t2nllT/N \m,l,s) , (7) 



with s being the index of two in-cell atoms and 2Mh the 
tubule length. In Eq. local state |m, I, s) of unit cell 
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(m, I) can be obtained from the (0, 0) unit cell |0, 0, s) by 
m successive screw operations S(a, h) combined with I 
successive rotations Cm as |m, I) = T m >|0, 0) with T m > = 
S m (a,h)C l N . 

For purpose of ORP calculation, we start from the off- 
diagonal element of the dielectric tensor s xy (q,u>), 



as A g = \/3ttVo/A. This gap sets a lower bound of the 
incident light wavelength as A c = 2hcA / (\/3Po) , below 
which no optical transverse absorption take place. Since 
this threshold is independent of chirality, it provides a 
means to determine the tubule diameter for all kinds of 
SWCNs by optical measurements. 



87re x -> x -> 
m* 2 oj 2 V ^ j-f 

Ktia K.'n'a 



(f(EZ,( K '))-f(EZ(K))) 



(nna\Q x (—q)\ K'n'a') (K,'n'o~'\Q y (q) \ una) 
hu} + EZ(K)-E%(K')+ie 



, (8) 



with a the band index c or v. Here, m* and e are the elec- 
tron mass and charge respectively. V = p ■ Mh\R\ 2 / (2tt) 
is the system volume with p as the filling factor. At zero 
temperature f(E) is simply the step function and the 
vector operator Q(g) has the form [p + eA/c, e tqz ] + /2 
whose x and y components transform under symmetry 
operation T m i as 




T m,l[Qx(l) ± iQy(q)]T m .i 
= e mqh±l( - ma+l7 '^[Q x (q)±iQ y (q)} . 

One then obtain the matrix elements of Q x ±iQ y between 
the conducting and valence bands to be 



(«' ',n' ' ,c\Q x (q) ±iQ v (q)\n,n,v) = 5 K ^ K+qh ± a 8 n ^ n ±i 
J^exp 



Y I . i \ . 2(n±l)7r , 
-i[k ± a + qhjm — i — i 



N 



x (m, I, c\Q x (q) ± iQ y (q)\0, 0, v) , 



(9) 



here the summation is over the nearest neighbors of the 
(0, 0) unit cell and it reflects the local geometric struc- 
ture of SWCN. The 5-function in Eq. © gives the se- 
lection rule for the transverse optical transition. To 
keep the gauge invariance, we transfer the calculation 
of the momentum matrix element into that of the co- 
ordinate matrix element through the commutation law 
p + eA/c = m*[r, H]/(i%)2Q to avoid the explicit treat- 
ment of vector potential A. The translational invari- 
ance is also rigorously kept in our calculating procedure. 
Based upon all the above considerations, the ORP can 
be calculated as 



X 



e 2 , , 2 v 2 
6 U) Vq \ 

n=l 



dts 



W(k, n)A(n, n)d K D(n, n) 



[(J%u> + ie) 2 — A(k, n) 2 ] 2 
7n(«0l£+i(« + ") - ivA 



TTpC 

W(k, n) = 1 - Re . , 

.|7n(«)||7n+x(« + Q!)| 
D(K,n) = |7„+i(K + a)| 2 - |7„(k)| 2 , 

A(«,n) = £?S 1 (« + a)-4^(«), 



(10) 



where tfAB = 7r(ni + ?i2)A _2 is the difference between 
the azimuth angles of the two atoms in one unit cell. 

In the optical limit q — > the minimum of A(k, n) gives 
an indirect band gap for the transverse optical transition 




FIG. 2: The dependence of the ORP on magnetic flux and 
chiral angle is plotted in (a) for v = 1 and in (b) for v = 0. 
The normalized incident light frequency is ui = 0.5 and the 
tubule diameter is around 37 nm. In (a) the ORP for u = 1 
tubules is a periodic function of chiral angle with period 2tt/3 
and magnetic flux with period <f>o- The ORP for v = — 1 
tubules can be obtained through the translation of chiral angle 
by 7r/3 in (a). In (b), the ORP for the v = tubules is a 
periodic function of chiral angle with period 7r/3 and magnetic 
flux with period <j)Q. 

When plotting \ versus the chiral angle 8 for arbitrary 
[77,1,712], we found the numerical data can be interest- 
ingly regrouped into three categories in accordance with 
the chiral index v — 0,±1. By introducing a rescaled 
dimensionless frequency Q — Tko/A g , x" can be fit- 
ted perfectly by the following functions with parameter 
/ = 4iu 2 e 2 /npc 2 Vo 



X ± /f « af(^,^-)sin(30)+A- 1 a 2 ± (w,^)sin(60) 
<Pa ' <Po 



X°/f 



A- 1 a^(w,^-)sin(60) , 
00 



(11) 



where af — —af and — a^ which, as dimensionless 
functions of Q and <p/(f>o, can be determined numerically. 
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Eq. Qllj l is entirely consistent with the above symmetry 
analysis (Eqs. J2J and J3J). In Fig. 2 we give 3D plots of 
X + and \° versus magnetic flux and chiral angle. Note 
that x~ can be obtained through exact relation x~(0) = 
X + {6 +7r/3). The results explicitly shows AB oscillation 
in the ORP with a single flux quantum period and x v is 
an even function of <f>. For large A, the numerical data 
verifies x + = ~X~ within its precision and x° is much 
smaller than x^ m magnitude (see also Eq. (111(1 ). It is 
interesting to note that for v = ±1 tubules the magnitude 
of x^ reaches its maximum/minimum when chiral angle 
approach 7r/6 for fixed <fi and u>, however for the armchair 
tubules (6 = 7r/6) x° — 0. This implies that the chiral 
index plays a role much more sensitive than the chiral 
angle in the natural gyrotropy properties of SWCN. 
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FIG. 3: Coefficients a±, aj, and a% are plotted as functions of 
magnetic flux <j> for different values of renormalized frequency 
uj in (a), (b), and (c), respectively. The filled square, open 
circle and cross symbol correspond to u — 0.4, 0.5 and 0.6, 
respectively. Fig.(d) is the plot of a+ , aj, and a% vs. Co de- 
noted by triangle, circle and square, respectively, for different 
values of magnetic flux. The open symbols are the data for 
cj> — and the filled symbols for cj> — 0.5(po- 

The fitting coefficients , aj, and are plotted in 



Figs. 3a, 3b, and 3c respectively as the functions of </>/</>o 
with Q = 0.4, 0.5, 0.6. All af with different Q coincide 
approximately at a common node 4>/4>o — 0.31, i.e., co- 
incide at ±0.31 + fi with fj, being the integer, since it 
is a periodic even function of 4>/4>q. However, for the 
v = tubules, the position of the zero point changes dis- 
tinctively versus uj. Since dominates the x function, 
this remarkable fact provides a way to distinguish exper- 
imentally whether the SWCN under investigation belong 
to v = ±1 or v = 0. In Fig. 3d we plotted , aj, a" as 
functions of u> for fixed <p/<t>Q = 0, 0.5. It is seen that when 
the incident light frequency approaches the band edge the 
ORP for both two kinds (u = and v = ±1 ) of SWCN 
increases rapidly. This fact may also be helpful to real- 
ize the ORP with enough intensity by tuning the proper 
frequency of the incident light. To have a quantitative 
estimation for the ORP of SWCN, we choose tubule di- 
ameter 10 A and the incident light wavelength 1.4 /im, 
= 0, and p = 2, which corresponds to Q ~ 0.8 and 
A ~ 13, then the magnitude of the rotation angle per 
unit length can be calculated to be 79 rad/cm x vsm(38) 
for v — ±1 tubules and 23 rad/cm x sin(60) for v = 
tubules. 



As the final remarks, the AB effect is conventionally 
investigated in connection with the transport studies for 
the cylindrical metallic sheet. In this paper, we show that 
the AB oscillation will also appear in the optical phenom- 
ena generically for arbitrary chiral SWCNs. Our argu- 
ments apply even to other kinds of chiral single wall nan- 
otubes. Moreover, the response to the enclosed magnetic 
flux is dramatically affected by the chiral index which 
characterizes different kinds of global helicity states of 
SWCNs and permits distinguished chirality dependence 
of physical properties. The gauge invariant TBA cal- 
culation not only provides quantitative results for var- 
ious quantities to stimulate the experimental measure- 
ments but also verifies all the conclusions drawn from 
our generic symmetry analysis. 
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